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Abstract Results from 2.5-D electromagnetic hybrid simulations show the formation of field-aligned, 
filamentary plasma structures in the magnetosheath. They begin at the quasi-parallel bow shock and 
extend far into the magnetosheath. These structures exhibit anticorrelated, spatial oscillations in plasma 
density and ion temperature. Closer to the bow shock, magnetic field variations associated with density and 
temperature oscillations may also be present. Magnetosheath filamentary structures (MFS) form primarily in 
the quasi-parallel sheath; however, they may extend to the quasi-perpendicular magnetosheath. They occur 
over a wide range of solar wind Alfvenic Mach numbers and interplanetary magnetic field directions. At lower 
Mach numbers with lower levels of magnetosheath turbulence, MFS remain highly coherent over large distances. 
At higher Mach numbers, magnetosheath turbulence decreases the level of coherence. Magnetosheath 
filamentary structures result from localized ion acceleration at the quasi-parallel bow shock and the 
injection of energetic ions into the magnetosheath. The localized nature of ion acceleration is tied to the 
generation of fast magnetosonic waves at and upstream of the quasi-parallel shock. The increased pressure 
in flux tubes containing the shock accelerated ions results in the depletion of the thermal plasma in these 
flux tubes and the enhancement of density in flux tubes void of energetic ions. This results in the observed 
anticorrelation between ion temperature and plasma density. 


1. Introduction 

The properties of the magnetosheath are determined by the conservation of mass, momentum, and energy 
across the bow shock (the Rankine-Flugoniot conditions) and collisionless dissipation processes upstream at 
and downstream of the bow shock. The most important solar wind parameters that impact the strength and 
the geometry of the bow shock are the Mach number, plasma beta (ratio of thermal to magnetic pressure), 
and the direction of the interplanetary magnetic field (IMF). The direction of the upstream magnetic field 
determines the regions of the bow shock corresponding to the quasi-parallel and quasi-perpendicular 
geometries, where shock dissipation processes are considerably different. In the quasi-perpendicular geometry, 
shock dissipation results in ions downstream from the bow shock having temperatures perpendicular to the 
magnetic field larger than parallel [e.g., Leroy et al., 1 982; Sckopke et at., 1 983, 1 990]. The presence of this 
anisotropy results in the generation of ion cyclotron and mirror mode waves in the magnetosheath with the 
former propagating primarily along the magnetic field and latter simply being carried by the magnetosheath 
flow, because they have zero phase and group velocities [e.g., Luhmann et al., 1 986; Lacombe et al., 1 992; 
Belmont et al., 1 992; Gary, 1 992; Gary et al., 1 993; McKean et al., 1 992, 1 994; Fazakerley and Southwood, 1 994; 
Anderson et al., 1 994; Omidi et al., 1 994; Lacombe and Belmont, 1 995; Lucek et al., 1 999, 2001 ; Genot et al., 2009]. 

In the quasi-parallel regime, ion reflection and leakage at the shock result in the presence of a variety of 
backstreaming foreshock ion beam populations [Asbridge et al., 1968; Greenstadt et al., 1968, 1980; Gosling 
et al., 1 978; Paschmann et al., 1 979; Bonifazi et al., 1 980a, 1 980b]. The interaction of these ions with the solar 
wind generates a number of different ULF waves and nonlinear structures that get carried back by the solar 
wind and impact the magnetosheath [e.g., Fairfield, 1969; Floppe et al., 1981; Russell and Floppe, 1983; Gary 
et al., 1 984; Le and Russell, 1 992; Greenstadt et al., 1 995; Blanco-Cano and Schwartz, 1 995; Lin, 2003; Omidi, 
2007; Blanco-Cano et al., 2009, 2011; Kajdic et al., 201 0, 201 1 ; Omidi et al., 201 3; Zhang et al., 2013]. The 
presence of tangential or rotational discontinuities in the solar wind and their interaction with the foreshock 
or bow shock result in the formation of structures such as foreshock bubbles or hot flow anomalies with 
significant effects on the magnetosheath [e.g., Schwartz et al., 1988, 2000; Schwartz, 1995; Thomsen et al., 

1 986, 1 988, 1 993; Paschmann et al., 1 988; Thomas et al., 1991; Sibeck et al., 1 998, 1 999, 2000; Lin, 1 997, 2002; 
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Facsko et al., 2008; Eastwood et at., 2008; 
Jacobsen et at., 2009; Omidi et at., 2010; 
Turner et at., 2013]. Solar wind transients 
and their interactions with the bow shock 
can also result in the generation of waves 
and density pulses in the magnetosheath 
[e.g., Sibeck and Costing, 1 996; Nemecek et at., 
1 998; Seon et at., 1 999; Zastenker et at., 2002; 
Savin et at., 2008; Amata et at., 201 1; Turner 
et at., 201 1; Karlsson et at., 2012]. 

In this paper, we use the results from 2.5-D 
global hybrid simulations to demonstrate for 
the first time the formation of field-aligned 
density and temperature enhancements that 
we call magnetosheath filamentary 
structures (MFS). These are large structures 
that begin at the quasi-parallel bow shock 
and extend well into the magnetosheath as 
far as the magnetopause. MFS form over a 
wide range of solar wind Mach numbers and 
IMF cone angles (defined as the angle 
between IMF and solar wind velocity). We 
show that their formation is tied to localized 
ion acceleration processes at the quasi-parallel 
bow shock and is a direct result of ion 
dissipation processes at the bow shock. In the 
following section, we briefly describe the 
hybrid simulation model and show the results 
that illustrate the formation of MFS under a 
variety of Mach numbers and IMF cone angles. 
Section 3 provides a summary and conclusions. 

2. Simulation Results 

We use 2.5-D (3-D in velocity and fields and 2- 
D in space) electromagnetic hybrid (kinetic 
ions and fluid electrons) simulations to 
examine the formation of MFS. This model 
has been used extensively in the past and 
has resulted in a number of predictions 
subsequently observed by spacecraft [e.g., 
Omidi et at., 2004, 2005, 2006, 2009a, 
2009b, 2010, 2013; Omidi and Sibeck, 2007; 
Blanco-Cano et al., 2006a, 2006b, 2009, 
2011; Sibeck et al., 2008; Turner et al., 2013; 
Zhang et al., 2013]. The four equations that are solved include Ampere's law without the displacement 
current (called the Darwin approximation) so that high-frequency light waves are eliminated. The 
second equation is Faraday's law, while the third is the electron momentum equation with electron mass 
equal to zero. This leads to the generalized Ohm's law that includes the Hall term and electron pressure 
gradients. Electrons are assumed to behave adiabatically. The model has a magnetic dipole inside a 
conducting sphere whose surface represents the ionospheric boundary. The simulation box lies in the 
noon-midnight meridian plane with x along the solar wind velocity (Sun-Earth line) and the magnetic 
dipole moment in the z direction. The simulation domain covers 2000 ion skin depths c/m p (where c is 



Figure 1 . Results from a run with Alfven Mach number of 3 and IMF 
cone angle of 35°. (a) The plasma density normalized to solar wind 
value and magnetic field lines, (b) The ion temperature normalized 
to its solar wind value, (c) The total magnetic field strength normalized 
to IMF value with the maximum value artificially set at 10 to make the 
magnetosheath more visible. 
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Figure 2. (a) The plasma density zoomed near the quasi-parallel bow shock and magnetosheath, (b) The ion temperature 
in the same area as in Figure 2a. Figures 2a and 2b both show density and temperature enhancements associated with MFS. 


the speed of light and co p is the ion plasma frequency) in the x and z directions with a cell size of 1 ion 
skin depth (see Omidi and Sibeck [2007] for the discussion on the choice of cell size in global hybrid 
simulations). Solar wind plasma is injected from the left boundary and allowed to leave the system 
through the remaining boundaries. The results discussed here correspond to solar wind conditions with 
Alfven Mach number between 3 and 11 and IMF cone angles (with respect to x axis) of 35° and 5°. 

Figure 1 shows the output from a run with an Alfven Mach number of 3 and IMF cone angle of 35° zoomed 
mostly on the quasi-parallel (southern) portion of the bow shock. Figure la shows the density normalized to its 
solar wind value and magnetic field lines. Figure 1 b shows the ion temperature normalized to the corresponding 
solar wind value, and Figure 1 c shows the total magnetic field strength normalized to the IMF value. In Figure 1 c, 
the scale is set to a maximum of 1 0 in order to make the bow shock and magnetosheath visible. The density and 
temperature panels show the presence of field-aligned, filamentary structures in the quasi-parallel 
magnetosheath. On the other hand, there are no signs of filamentary magnetic field structures within the 
southern magnetosheath behind the quasi-parallel bow shock where, except near the shock, magnetic field 
strengths are comparable to those in the solar wind as expected from the Rankine-Hugoniot conditions. 

To view the magnetosheath filamentary structures in more detail, Figure 2 shows the density and temperature 
from the same run as in Figure 1 in the quasi-parallel portion of the sheath. It is evident in Figure 2b that the 
filamentary structures in the magnetosheath are connected to the regions of localized high-energy ions 
upstream of the shock. Figure 3 shows the variations of the magnetic field, ion temperature, and density along 
the lines LI and L2 (length of ~1 00 c/co p ) in Figure 2a. Figures 3a-3c correspond to LI and show an evidence of 
oscillations in the magnetic field, temperature, and density. It is evident that the wavelength of the oscillations 
in temperature and density is shorter than that associated with the magnetic field. The latter is associated with 
waves in the magnetosheath not connected to MFS. In addition, the fluctuations in density and temperature are 
anticorrelated. Figures 3d-3f correspond to line L2 and again show anticorrelated fluctuations in density and 
temperature with no corresponding variations in the magnetic field. The results in Figure 3 indicate that the 
magnetosheath filamentary structures are primarily a nonmagnetic phenomenon; however, as we show below 
near the bow shock, these structures may be associated with magnetic field oscillations. 

Figures 4a through 4f show the ion temperature at 6 times Of = 21 2.5, 275, 337.5, 400, 462.5, and 500 
(where L2 is the proton cyclotron frequency) and illustrate the evolution of the magnetosheath filamentary 
structures. MFS initially form at lower latitudes closer to the nose of the bow shock and subsequently get 
carried away from the nose as they continue to evolve. As a result, the number of filamentary structures 
increases with time as new ones are generated near the nose and then converted antisunward to higher 
(southern) latitudes. Figure 5 shows the transverse or dawn/dusk (B y ) component of the magnetic field at 
the same times as in Figure 4 and demonstrates the presence of waves upstream of the bow shock at all 
times shown. These are circularly polarized waves propagating along the magnetic field and are generated 
by the backstreaming ions seen in Figure 4 (see Blanco-Cano et al. [2006a, 2006b, 2009, 2011], Omidi 
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Along Line LI Along Line L2 

Figure 3. (a-c) The variations of magnetic field strength, ion temperature, and density along line labeled LI in Figure 2a. 
(d-f) The variations of magnetic field strength, ion temperature, and density along line labeled L2 in Figure 2a. Note the 
anticorrelated oscillations in density and temperature associated with MFS. 
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Figure 4. (a-f) The ion temperature in the foreshock and quasi-parallel magnetosheath at different times illustrating the 
formation and evolution of MFS. 
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Figure 5. (a-b) The transverse component of the magnetic field at the same time as those in Figure 4. They show the pre- 
sence of waves in the foreshock at all times and also the generation of oscillations in the magnetosheath near the quasi-parallel 
bow shock. 


[2007], and Omidi et at. [2005] for detailed analysis and discussion of upstream waves in global hybrid 
simulations). They have wave vectors that point toward the Sun, but because their group velocity is 
smaller than the solar wind speed, they are convected back toward the shock. Figure 5 also shows the 
presence of waves in the magnetosheath near the bow shock that, similar to MFS, originate close to the 
nose of the bow shock and travel to higher latitudes with time. The connection between MFS and these 
waves will be discussed shortly. 



Figure 6. (a-b) The total magnetic field at the same time as those in Figures 4 and 5. They show the generation of fast mag- 
netosonic waves at large wave angles with respect to magnetic field. 
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Figure 7. (a-d) The density, ion temperature, transverse magnetic field, and total magnetic field zoomed near the quasi-parallel 
bow shock and sheath. Also shown are the same set of magnetic field lines in each panel. 




Along Line L3 Along Line L4 

Figure 8. (a-d) The variations of ion temperature, density, magnetic field strength, and transverse magnetic field along line 
labeled L3 in Figure 7c. (e-h) The variations of ion temperature, density, magnetic field strength, and transverse magnetic 
field along line labeled L4 in Figure 7c. 
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Figure 9. (a and b) The plasma density and ion temperature 
zoomed near the quasi-parallel bow shock from a run with Alfven 
Mach number of 5 and IMF cone angle of 35°. 


Figures 6a through 6f show the total 
magnetic field strength at the same times as 
in Figures 4 and 5. Early in the run (Figure 6a), 
no fluctuations in the total magnetic field are 
present in the foreshock, consistent with the 
transverse nature of the parallel propagating 
waves in that region. At later times, Figure 6 
shows the appearance of magnetic 
fluctuations near and upstream of the bow 
shock that originate at lower latitudes and 
appear at higher latitudes with time. These 
are fast magnetosonic waves that propagate 
at large angles to the upstream magnetic 
field. They were first discussed by Omidi 
[2007] and Blanco-Cano et al. [2011], who 
showed that the backstreaming ions in the 
foreshock not only generate the transverse, 
parallel waves but also highly oblique fast 
magnetosonic waves. The growth rate of 
these waves increases with the energy of the 
backstreaming ions. So in the initial formation 
of the shock in the simulation, the energy of 
the backstreaming ions is not high enough for 
appreciable growth of magnetosonic waves, 
and only parallel propagating waves are 
excited (Figure 6a). Flowever, as higher-energy 
ions originating from the nose region of the 
magnetosheath leak into the upstream 
region, fast magnetosonic waves are also 
excited (Figures 6b-6f). These waves are 
converted into the magnetosheath and result 
in the oscillations in B y seen in Figure 5. 


Magnetosonic waves are responsible for the localized acceleration of the ions at the bow shock and the 
formation of MFS. Figure 7 illustrates more details on the structure of MFS and shows density, temperature, 
B y , and total magnetic field zoomed on the quasi-parallel bow shock and magnetosheath. Also plotted are 
magnetic field lines in the four panels so that we can focus on individual field lines and examine the details of 
the MFS structure. For example, the examination of the third magnetic field line (thicker and green color) 
from the upper left corner in Figure 7a shows that it falls in a region of low density in the magnetosheath. 
Figure 7b shows that the same field line falls in a region of high ion temperature in the magnetosheath 
consistent with the anticorrelation of density and temperature described earlier. It is also evident that this 
field line connects to a region of higher ion temperature in the foreshock and at the bow shock. As can be 
seen in Figure 7d, this region of higher temperature corresponds to lower magnetic field strengths and is 
bounded by higher magnetic fields associated with the fast magnetosonic waves. Similar patterns are 
associated with the other field lines depicted in Figure 7, indicating that the magnetosonic waves cause 
localized ion temperature enhancements in the foreshock. The energized ions follow magnetic field lines into 
the magnetosheath. Flere enhanced pressures associated with the energetic ions cause expansions that 
diminish plasma densities on these field lines. By contrast, neighboring regions with fewer or no energetic ions 
are compressed, resulting in density enhancements. This leads to the consistent anticorrelation seen between 
density and temperature changes associated with MFS. 


The examination of Figures 7c and 7d shows that the fluctuations in B y and total magnetic field are present in 
the sheath near the bow shock. The regions of enhanced B y and field strength coincide with each other and 
also line up with peaks in the magnetic field strength upstream of the shock. This indicates that the 
magnetosheath oscillations result from the convection of upstream magnetosonic waves into the 
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downstream. Since the structure of the MFS 
and the spacing between the filaments are 
controlled by the upstream magnetosonic 
waves, one would expect some form of 
correlation between the oscillations in B y and 
total field and MFS near the bow shock. To 
examine this possibility, Figures 8a-8d show 
ion temperature, density, magnetic field 
strength, and B y along the line L3 shown in 
Figure 7. It is evident that density and 
magnetic field are correlated with each other 
and anticorrelated with temperature. The 
correlation between density and magnetic 
field is consistent with the signatures 
expected for a fast magnetosonic mode. 
Figures 8e-8h show the variations of the 
magnetosheath plasma and magnetic field 
along the line L4 in Figure 7 and demonstrate 
that further away from the bow shock, MFS 
shows no correlated changes in magnetic 
field. Thus, given the connection between the 
upstream fast magnetosonic waves and MFS 
and the convection of the former into the 
magnetosheath, we find oscillations in the 
magnetic field that correlate with MFS near 
the bow shock. Flowever, further into the 
magnetosheath, the MFS has no magnetic 
component and consists solely of density and 
temperature oscillations. 

The results of hybrid simulations show that 
MFS form over a wide range of Mach 
numbers and IMF cone angles. Figures 9a and 
9b show the density and temperature from a run with IMF cone angle of 35° and solar wind Alfven Mach 
number of 5. The formation of MFS in the quasi-parallel magnetosheath is evident and demonstrates that it is 
not only associated with low Mach number bow shock. There is however a change in the quasi-parallel bow 
shock at higher Mach numbers in that spontaneous hot flow anomalies (SFHFAs) are formed as part of the 
shock dissipation process [Omidi et at., 201 3; Zhang et at., 201 3]. SFIFAs are cavities containing low solar wind 
densities and magnetic field strengths and enhanced fluxes of high-energy particles with rims exhibiting 
enhanced densities and magnetic field strengths. They form as a result of the convection of foreshock 
cavitons into the quasi-parallel bow shock. The formation of foreshock cavitons was predicted by hybrid 
simulations [Lin, 2003; Lin and Wang, 2005 ; Omidi, 2007], and they were discovered in spacecraft data by 
Blanco-Cano et ai. [2009, 2011] and Kajdic et at. [2010, 2011]. As demonstrated by Omidi [2007], foreshock 
cavitons result from the simultaneous generation of parallel propagating electromagnetic waves and oblique 
fast magnetosonic waves by the backstreaming ions and their highly coupled nonlinear evolution. The 
formation of SFIFAs at the bow shock also results in localized regions of high-energy ions that are injected 
into the magnetosheath along the field line resulting in the formation of MFS. We note that at both low and 
high Mach numbers, the role of fast magnetosonic waves in forming localized regions of higher ion 
temperature at the bow shock is similar. 

The results presented so far correspond to an IMF cone angle of 35° with MFS forming mostly away from the 
dayside magnetopause and carried away with the magnetosheath flow. At lower IMF cone angles, the MFS 
forms upstream of the dayside magnetopause and stays more stationary due to the low flow velocity of the 
magnetosheath in that region. Figures 10a and 10b show the density and ion temperature from a run with 
an Alfven Mach number of 5 and an IMF cone angle of 5°. The formation of MFS in the low-latitude dayside 
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Figure 10. (a and b) The plasma density and ion temperature from a 
run with Alfven Mach number of 5 and IMF cone angle of 5°. 
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magnetosheath is evident in both the 
density and temperature plots. In contrast 
to the case for higher IMF cone angles, 
where the filaments lie nearly parallel to 
vector pointing radially away from Earth 
and move poleward, the density and 
temperature perturbations associated with 
the MFS in Figure 1 0 lie nearly perpendicular 
to the radial direction. Figure 11 is similar 
to Figure 10 except that it corresponds to a 
run with a higher Alfvenic Mach number of 
11. The MFS structure in Figure 11 is similar 
to that in Figure 10, because the IMF cone 
angle is the same in both runs. Comparing 
Figure 10a and Figure 11a demonstrates 
that while MFS forms in both runs, the 
magnetosheath in Figure 1 1 has a more 
turbulent structure due to the higher solar 
wind Mach number. The enhanced 
turbulence causes the magnetic field lines 
that comprise the MFS to bend and twist 
into more complicated forms. 

3. Summary and Conclusions 


Figure 1 1 . (a and b) The plasma density and ion temperature from a 
run with Alfven Mach number of 1 1 and IMF cone angle of 5°. 


Using global, electromagnetic hybrid 
simulations, we show the formation of 
magnetosheath filamentary structures: 
namely, field-aligned density and 
temperature structures in the 
magnetosheath. They begin at the quasi- 
parallel bow shock and extend well into 
the quasi-parallel magnetosheath reaching as far as the magnetopause. When the IMF cone angle is 
small, MFS may extend into the quasi-perpendicular magnetosheath. The antisunward convection of 
MFS that are stationary in the plasma rest frame results in anticorrelated magnetosheath density and 
temperature fluctuations. The fluctuations are nearly stationary in the low flow speed subsolar 
magnetosheath but produce Doppler shifted frequencies in the higher-latitude magnetosheath, where 
flow speeds are greater. In general, MFS are not associated with magnetic field changes; however, near the 
bow shock, the presence of transmitted fast magnetosonic waves may show magnetic field oscillations that 
correlate well with MFS associated density and temperature variations. The convection of the parallel 
propagating waves in the upstream results in magnetic field variations at and downstream of the shock 
that would not correlate with density and temperature oscillations associated with MFS. 


Magnetosheath filamentary structures form over a wide range of solar wind Mach numbers and IMF 
cone angles. In all cases, their formation is tied to the presence of localized regions of enhanced ion 
temperature at and upstream of the quasi-parallel bow shock. Upon injection of these energetic ions 
into the magnetosheath, they follow the field line and form a temperature filamentary structure whose 
enhanced pressure reduces the plasma density in the same flux tube but increases it in the surrounding 
flux tubes. This accounts for the observed anticorrelations between density and temperature 
fluctuations associated with MFS. The formation of localized regions of enhanced ion temperature at 
and upstream of the bow shock is connected to the generation of fast magnetosonic waves propagating 
at large angles to the upstream magnetic field. At low Mach numbers, these waves result in local 
acceleration of ions upstream of the bow shock through interaction with the parallel propagating, 
circularly polarized electromagnetic waves that are also generated by the back streaming ions. At higher 
Mach numbers, the nonlinear evolution of the fast and parallel propagating waves results in the 
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formation of foreshock cavitons, whose convection into the bow shock results in the generation of 
SHFAs and localized enhancements of temperature. 

The scale length associated with MFS, namely, the distance between the filaments is primarily determined by 
the distance between the regions of high ion temperature at the quasi-parallel bow shock. This spacing is 
related to the wavelength of the fast magnetosonic waves and is -30-60 ion skin depth, which corresponds 
to -0.5-1 R e . For magnetosheath velocities of 1 00-300 km/s, the expected frequency of the MFS oscillations 
would be in the range of 0.01-0.1 Hz assuming flow along the wave vector. In general, this assumption 
may not be true, and a determination of the local angle between the flow velocity and the direction of the 
wave vector may be needed in order to assess the Doppler frequency of MFS oscillations. 
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